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A coherent study of the production of /q (i = 1, 2, 3 corresponding to /o(1710), /o(1500), and 
/o(1370)) in J/ip — > Vfo — > VPP is reported based on a previously proposed glueball and QQ nonet 
mixing scheme, and a factorization for the decay of J/tp — > Vfo, where V denotes the isoscalar vector 
mesons <p> and u>, and P denotes pseudoscalar mesons. The results show that the J/ip decays are very 
sensitive to the structure of those scalar mesons, and suggest a glueball in the 1.5—1.7 GeV region, 
in line with Lattice QCD. The presence of significant glueball mixings in the scalar wavefunctions 
produces peculiar patterns in the branching ratios for J/ip — > Vfo — ► VPP, which are in good 
agreement with the recently published experimental data from the BES collaboration. 



I. INTRODUCTION 



Although explicit evidence for a pure glueball state has never been confirmed in the spectroscopy of isoscalar 
mesons, experimental data are consistent with the scenario that a glueball should have significant mixings with 
quarkonia in the formation of some of those isoscalar mesons. This is expected to be especially important in the scalar 
mesons . It is hence interesting to explore those mechanisms of isoscalar meson production which are sensitive to 
their glueball contents. Recently, Brodsky et al. 2} proposed a possible production mechanism for glueball states 
in e + e~ — * 7* — > J/ijjQj, where glueball production in the same mass region as the scalar charmonia {n c and Xco) 
would contribute to the events observed in e + e~ — * 7* — * J/ipX at Belle |3J. In contrast to the above prediction for a 
glueball of about 3 GeV, Lattice QCD 00 and other phenomenological studies 0,0,13 suggest the lightest glueball 
mass at 1.5 ~ 1.7 GeV. Some analyses identify the lowest + glueball with the /o(1500) 0,13,0], while even lower 
glueball masses have been suggested in the literature [ll|. It is thus natural to concentrate on this energy region 
and explore glueball production and possible mixings with quarkonia in the low-lying scalars, in particular, the three 
states /q, where i = 1, 2 and 3 correspond to /o(1710), /o(1500) and /o(1370), respectively. 

In Refs. , Close et al. proposed a flavor mixing scheme to produce the three /o states among the glueball G, 
flavor nonet ss and nn. It hence relates the decays /q — » PP, where P denotes pseudoscalar mesons (it, K, 77, 77'), to 
the mixing angles between G, ss and nn. Since the mixing scheme is sensitive to the experimental data for /q — > PP, 
their results not only provided a prescription for understanding the data from WA102 0], but also highlighted the 
importance of the possible glueball contents in both /o(1500) and /o(1710), while the /o(1370) was dominated by the 
nn component. In this scheme, the pattern of the mass matrix will be determined by the relative mass positions of 
G, ss and nn. Therefore, for the pure glueball G it is crucial to compare with the Lattice QCD predictions. 

Independent production processes are required to examine such a glueball-QQ mixing mechanism and provide more 
decisive evidence for the existence of the glueball contents in some of the /o states. At the Beijing Spectrometer 
(BES) the /o states have now been observed in both J/ip radiative and hadronic decays, and interesting signals of 
these scalars have been reported 0,0,0]. In particular, the data show that these states behave contrary to ex- 
pectations based on naive application of the Okubo-Zweig -Iizuka (OZI) rule in J/ip decays fift: i) the / (1370) 
has been seen clearly in J/tp — > (pTnr, but not in J/ip — > unrTt; ii) there is no peak of the /o(1500) directly seen in 
J/tp — > 4>KK, ujKK, 07T7T, ujirn, though its production in proton-proton scattering is quite clear 0; iii) the /o(1710) 
is observed clearly in both J/ip -> cj)KK and J/ip -> uKK, but with br J/4 ,^ ujMl7w) ^ ujKR /brj / ^^ ct>fo{1710) ^ 4lKIt ~ 6, 
which is against a simple ss configuration for this state. The /o(1710) is not seen in <^7T7r and lvttt:, which is under- 
standable due to its small branching ratio to tttt. These observations hence raise questions about the configuration of 
the /o states, and also about their production mechanism in the J/tp hadronic decays. BES also claims observations 
of another scalar /o(1790), which is distinct from the /o(1710). We will comment on this state in the last Section, 
based on our results for the /q' 2 ' 3 - 

An estimate of the low-lying glue ball production rate has been made by extrapolating the pQCD approach of 
Ref. to the J/ip energy region |l8j . By studying e + e~ — > J/tp — > <pG there, we showed that the glueball production 
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rate can be approximately normalized in terms of the ideal ss(lS) state production. One of the results of the present 
paper will be to extract from the data the effective branching ratios for J/ip — > VG. 

In what follows, first we shall revisit the glueball-QQ mixings in these three scalar /q mesons, and clarify quantita- 
tively the sensitivity of the mass matrix pattern to the available experimental data from WA102 ^1 

and BES [3[l|. 

We will then explore the implications of the glueball-QQ mixings in the /g productions in J/ip — > V/q, and compare 
the results with the recent experimental data from BES 0, Conclusions about the /q properties and their 

"puzzling" behaviors in the J/ip decays will then be drawn. 



II. THE THEORETICAL MODEL 

A. Glueball and QQ mixing 

References proposed a flavour mixing scheme for the three scalars /q' 2 ' 3 corresponding to /o(1710), /o(1500) 

and /o(1370), respectively. In the basis of \G) — \gg), \ss) and \nn) = \uu + dd) / \/2, the mass matrix for the 
glueball-quarkonia mixing can be written as 
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where Mg, M sS and M nn represent the masses of the pure states |G), \ss), and \nn), respectively. The parameter 
/ is the flavour independent mixing strength between the glueball and quarkonia, and leads to the mixture of \G), 
\ss), and \nn) to form the physical /q states. It can be expressed as / = (sslV^G) = (nn\V\G) /y/2, where V is the 
potential causing such a mixing. 

Treating /q(1710), /o(1500) and /o(1370) as the eigenstates of the mass matrix M with the eigenvalues of the 
physical resonance masses Mi, M2 and M3, respectively, we can then express these fo states as 

|/o(1710)) \ fx iyiZl \f |G) \ 

|/o(1500)) = I x 2 y 2 z 2 ] [ \ss) , (2) 
|/o(1370)) J \ x 3 y 3 z 3 J \ \nn) J 

where Xi , yi and Z{ are the mixing matrix elements determined by the perturbative transitions 0, Q • 
The decay of /q — > PP can then be factorized into contributions from three fundamental processes 8] : 

i) the direct coupling of the quarkonia component of the three fo states to the final pseudoscalar meson pairs; 

ii) the coupling of the glueball component of the three fo states to the pseudoscalar meson pairs; 

iii) the direct coupling of the glueball component to the gluonic component of the final state isoscalar meson pairs. 
The coupling strength of (ii) is denoted by parameter r 2 relative to the coupling of (i); analogously, parameter 7*3 is 

introduced to account for the coupling strength of (iii) relative to (i). In constrast to Ref. [Sj, we make here the more 
natural assumption that the G component of fo states couples to the flavour singlet QQ component of the 7/ and 77'. 
Taking into account the 77-77' mixing, the factorization for /q — > PP is presented in Table [I] 



B. The production of /o states in J/ip — > Vfo 

By defining the transition amplitudes via a potential V^, 

M 0G = {G\Vi\J/il>) 
M^s-s) = (ss|V^|J/V>) 

M 0(nfi) = {nn\V^m) , (3) 

we can then express the transition amplitudes for J ftp — > (pf^ in terms of the production of glueball |G), nonet \ss) 
and \nn), via their configuration mixings. Thus, the transition amplitudes can be factorized: 

Mj,w = {fl\G){G\Vt\J/^) + (folssXsslVtm) + (/S|nn)<nn|V,|JM 

= XiM^c + yiM^ sS ) + ZiM^njQ . (4) 

In Fig. ^ the leading diagrams of these transitions are illustrated. For J/ip — > <j>fg, the subprocesses (a) and (b) 
are of the same nominal order in perturbative QCD. Thus, we assume (G\V^,\J/ip) ~ (ss\Vj,\J / ip) . The transition 



-^0(rin) = (nn\V,f,\J/ip) can only occur via the doubly OZI disconnected process (c), which will be of 0(a s ) relative to 
(a) and (b) in the perturbative regime. However, due to the complexity of the hadronization of the gluons into QQ, 
the role played by the doubly OZI disconnected diagrams is still not well-understood. In particular, the strong mixing 
of nh-ss in scalar ++ illustrated by the sizeable mixings (see later Section) shows that Fig. may be expected 
to occur at similar strength to the singly disconnected diagram, Fig. ^b). To estimate its contribution, we then 
introduce another parameter r, which describes the relative strength between the doubly OZI disconnected process 
and the singly disconnected one such that: 

(GTO J/ip) ~ (88\V^\ J/il>) a -^-(nn\V^m) , (5) 

where the factor y/2 is from the normalization of nh. In the above equation, we have eventually assumed that the 
potential is OZI-selecting. The left equation holds since (G\V^\J/ip) and (ss\V^\J/ip) both are of the same nominal 
order and both are singly OZI disconnected. The right equation then distinguishes the doubly OZI disconnected 
process from the singly disconnected one. It should be noted that contributions from the doubly disconnected processes 
with the same flavors, such as (pfo(ss) (e.g., replace the nh in Fig.^c) with ss), could also contribute in the transition 
(ss|V^| J/ip)- However, following the OZI-selection assumption for the transition potentials, such contributions can be 
regarded as being absorbed into the potential strengths. We will leave the r to be determined by experimental data. 
Naively, r << 1 would imply that the doubly disconnected process has a perturbative feature, while r ~ 1 could be 
a consistency check for the non-perturbative dominance of nh-ss mixing in the ++ channel |l9L |2CI| . 
The above assumption will then allow us to express the partial decay width as 

r j/4,^4,f^<t,pp = j^^br f ^ P p[x t +y t + V2rz i } 2 Tj / ^^ ct>G , (6) 

where p^i and p^G are the momenta of the <f> meson in J/ip — > 0/q and the virtual J/ip — ► (pG, respectively. Their ratio 
gives the kinematic correction for the decays of J/ip to the (p meson and states with different masses, i.e. Mi =^ Mq. 
Tj/^-^^g = Co|P0g||M/>g| 2 is the J/ip decay width into (p and glueball G, and \M^g\ 2 is the invariant amplitude 
squared; Co is the phase space factor. 

Similarly, the transition amplitudes for J/ip — > lo/q can be obtained: 

Mj/^fi = (r \G)(G\V u \J/1>) + (f \ss)(ss-K\JM + (mnn)(nn\VMm 

= XiMuG + ViM^gg) + ZiM^nn) ■ (7) 

Note that here the doubly OZI disconnected process is M w ( s5 ) = (ss\V u \J/ip). Also, for the exclusive process J/ip — > 
<x>/q — > luPP, we have the partial decay width 

r j/^ u /-^PP = r^br f ,^ P p[x i +ry l + V2z l } 2 Tj / ^ ujG , (8) 

where Tj/^^^g = CqIp^gI |Mjg| 2 is the J/ip decay width into uj and pure glueball G; p^^ and p^G represent the 
three momenta of the u> meson in J/ip — > w/q and J/ip — > loG, respectively. The flavour-blind assumption implies 
|-W wG | 2 = 2|M 0G | 2 , or equivalently, 2r J/v ,^0 G /|p0 G | = Vj/^ uG /\Pu)g\- For given /g, this relation will allow us to 
relate the oj and <p channel together, and then the parameter r can be determined. 



III. NUMERICAL RESULTS 
A. Parameters and the /q mass matrix 

In this Section, we will quantify these subprocesses and investigate the manifestation of the glucball-QQ mixings in 
the /q production mechanism in the J/ip decays. As summarized in the Introduction, the recent data for /q production 
in ip decays exhibit apparent differences from the simple QQ configurations. It is hence crucial for us to examine that 
the glueball-QQ mixings obtained in the study of /g — > PP are consistent with the observations in the J/ip decays. 
Predictions for the /g productions in J/ip — > V /g — ► VPP will then provide a self-consistent check of this approach. 

We first determine the model parameters by fitting the WA102 data Taking into account the unitarity and 
orthogonality constraint for the mixing matrix, the data can be fitted with a low reduced \ 2 ( see the Fit-I of Tabic [Hj, 
and the parameters are presented in Table ITTT1 Note that the /g decay branching ratio fractions generally have large 



errors. It shows that the fraction 6r(/ (1710) -> 7T7r)/6r(/ (1710) -> KK) = 0.20 ± 0.03 with the relatively small 
errors is the major contribution to % 2 . We also note that the constraints of unitarity and orthogonality will add 
another 9 relations in the fitting. We arbitarily require the error to be 3%. Because of this treatment, the absolute 
values of the x 2 only provide a reference for the fitting results. 

As shown by column Fit-I in Table ITT1 the masses of the |G), \ss), and \nn) follow the ordering of M sS > Mq > M nn , 
which is the same as Ref. 0, Furthermore, in this approach we find that the 77 — 77' mixing angle, <£> — —18.5 ± 3.1, 
fits naturally well the expectations from linear mass relations [2l| . It is interesting to see that the sign pattern of the 
wavefunction mixing matrix as observed in Ref. Q is reproduced: 

|/o(1710)) = 0.39|G) + 0.91|ss) + 0.13|rm) 
|/o(1500)) = -0.73|G) + 0.37|ss) - 0.57|nrl) 

|/o(1370)) = 0.56|G) -0.12|ss) - 0.82|rara) , (9) 

where |ss), |G), and \nn) are dominant components in /q' 2,3 , respectively. 

The BES experiment measured the exclusive branching ratios for the /q productions in four channels, J/iJj — > <j)KK 1 
(jmit, ujKK and lotttt. In particular, according to BES H|, 6r(/ (1710) -> 7T7r)/6r(/ (1710) -> KK) has an upper 
limit of 0.11 with 95% confidence level. This ratio, as pointed earlier, is an important contribution to the \ 2 m the 
fitting of the WA102 data. For the / (1370), the ratio 6r(/ (1370) -> K K)/br{f (1370) -> tttt) = 0.08 ± 0.08, is 
reported, which implies a much smaller KK decay rate than its tttt decays. In comparison with the WA102 data, 
these two ratios provide additional constraints. For the /o(1500), we derive &r(/o(1500) — > KK) /br(fo(1500) — * tttt) — 
(0.8 ± 0.5)/(1.7 ± 0.8) ~ 0.78 ± 0.66 by taking the bound values of the numerator and denominator. Although the 
value overlaps with the WA102 results (0.32 ± 0.07), the large errors will not be helpful for improving the \ 2 m the 
fitting. We hence include both in the fitting and the parameters and results are presented as Fit-II in Table [n] and 
Table ITTT1 to compare with Fit-I. 

Again, we find the consistent 77 — 77' mixing angle, <f> — — 22.7±3.5 |2l|. The same mass ordering, M sS > Mq > M nfi , 
and sign pattern of the mixing matrix are also preserved: 

|/o(1710)) = 0.36|G) +0.93|ss) + 0.09|nn) 
|/o(1500)) = -0.84|G) + 0.35|ss)-0.41|nn) 

|/o(1370)) = 0.40|G) -0.07|ss) -0.91|nn) , (10) 

which suggests the consistency between those two sets of data. 

To be specific, we abstract some stable features arising from these two fitting results as follows: 

i) Our analysis shows a strong preference of the 77-77' mixing angle to be around —23° ~ —18° degrees, which is 
consistent with previous data |21|, and should be a good reference for this factorization scheme. 

ii) In the proposed mixing scheme the glueball couplings into isoscalar meson pairs turn out to be strong. This 
process has a strength of r% ~ 1 relative to the QQ couplings to the pseudoscalar meson pairs. This suggests the 
glueball-QQ couplings in the scalar channel at these masses are not perturbative and cautions against a naive pQCD 
interpretation of the relative importance of Figs. D^a), (b) an d ( c ), a t least for ++ mesons. 

iii) The same mass ordering, i.e. M sg > Mq > M nfL , is found in both of these fitting results as shown in Table UTTI 
This is in agreement with the results of Ref. 0- The mass for the physical /o(1370) is treated as a free parameter, 
and we find in Fit-I that M 3 = 1275 ± 34 MeV and in Fit-II, 1260 ± 13. The preference of low masses for the / (1370) 
may be consistent with a recent analysis 0, |z| . The analyses of Ref. suggest that the difficulty of separating the 
/2(1270) peak from the /o(1370) led to the difficulty of determining the mass pole position of /o(1370). Due to the 
rapid rise of /o(1370) — ► Att phase space with mass, a peak occuring at about 1390 MeV in the Att channel is actually 
farther away from the mass pole compared with the tttt channel. Refined analyses show that a pole position around 
1315 MeV may be reasonable [22|. 

iv) In both fittings the masses of the pure \ss) are found to be lower than the physical /o(1710), which is dominated by 
the \ss) component. Furthermore, the pure \nn) are found to be higher than the physical /o(1370), which is dominated 
by the \nn). This suggests that the glueball mixings lead to a larger mass gap between /o(1370) and /o(1710) in 
comparison with the mass gap between the simple \nn) and \ss). This is also consistent with the observations that 
the mass gaps between the QQ dominated axial-vectors, /i(1285) and /i(1530), and between the tensors, ^(1270) 
and /2(1525), are relatively smaller than that between /o(1370) and /o(1710). Furthermore, the relative elevation of 
the mass of the /o(1710) relative to the /i,2(~ 1530) states, in contrast to the similar mass for the low mass states 
fo, 1,2(1270 ~ 1370) suggests that the primitive glueball lies in the 1.5 ~ 1.7 GeV mass range (this is also in agreement 
with our analysis of the scalar states, which follows). 

v) The most interesting feature is the stable pattern of the wavefunction mixing matrix arising from these two 
fittings, which is essentially determined by the relative mass positions between the pure states (|G), |ss), and \nn)) 



and physical ones (Mi. 2,3). Note that the mixing elements are given by 



Xi = (Mi- M sS )(M l - M nn )d 

Vi = (Mi — M nn )fCi (11) 
Zl = V2(Mi - M sM )fd , 

where Cj is the normalization factor of Eqs. © and l|l()|l. and both Cj and / are positive. For the /o(1370), as 
an example, the phases can be analytically understood with M sS > M nfi > M3. The same analysis applies to the 
/o(1710), where it is essential to have Mi > M sS . For the /o(1500), it shows that Mi can either larger or smaller than 
Mg if both are located between M nft and M sS . This also explains the relative strengths between those components 
in the wavefunctions. In particular, all the components are in phase for the /o(1710); the \ss) and \nn) are out of 
phase in the /o(1500) wavefunctions ("flavor octet tendency"); for the /o(1370), the \nn) and \ss) are in phase and 
hence have a "flavor singlet tendency". For the last, the significant \nn) and negligible \ss) turn out to agree with the 
observation that the /o(1370) has large branching ratios to irir and 4ir channels. 

In brief, we find that the inclusion of the recent BES data does not change the phase patterns of the wavefunction 
mixing matrix. The stable phase patterns and relative strengths between the glueball and QQ components can be 
understood as the preference of the experimental data to have M sS > Mg > M nfL with 1.46 < Mg < 1-52 GeV. 
This turns out to be meaningful since such stable mixing patterns possess predictive power, which can be tested in 
independent processes. 



B. Constraints from the /o(1710) production 

The mixing matrices provide a starting point for considering the exclusive decay processes of J/ip — * ^fo ~ * <fiPP 
and J/ijj — > lv/q — > ujPP. In particular, the clear signals for / (1710) in J/tp — + 0/ o (171O) — > 4>KK and J/ip — > 
o;/o(1710) — > ujKK channels 0,0] will allow us to normalize the glueball production rates via Eqs. © and (JSJ) in 
Section IIB. We will show in this Subsection that the /q production in J/tp decays is correlated with the glueball-QQ 
mixings introduced in the previous Subsection, and a self-consistent treatment of quark-glue couplings is needed. We 
then clarify that the BES data for the /o(1710) will provide not only a crucial constraint on the model, but also a 
self-consistent examination of this approach. 

From Eqs. 10 and (JSJ, we have 

Tj/y-^fi^KK = j^^j- ^-.KK^i + 2/1 + V2rzi] 2 r, 7/ ^_ >0G 

v j/i,^utt^uKK = \~\ br f^KR[xi +ryi +V2z 1 } 2 r J/i ,^ LoG ■ (12) 

IPwGl 

The BES collaboration 0| reported the partial decay branching ratio Tj/^^fi^xg/T^,^ — (2.0 ± 0.7) x 10~ 4 , 

which is consistent with the estimate of the PDG [2l| of (3.6 ± 0.6) x 10~ 4 . For T J/i> ^ ujf i^ KR /T T J/i: , BES [3 

found much larger ratio (13.2 ± 2.6) x 10~ 4 than the PDG [2l]], (4.8 ± 1.1) x 10~ 4 , due to better measurement of the 
J/ip — > cuKK channel. 

As discussed previously, the coupling r determines the role played by the OZI rules for the doubly disconnected 
processes. It is, in principle, independent of the /q decays, but should be determined by the /q production mechanisms 
in J/tp — > Vfy. To see more clearly the contribution of the doubly OZI disconnected subprocess, we take the 
naive OZI constraint, r — > 0, in the perturbative limit. Follow the implication of the flavour-blind assumption: 
2Fj/^,^0 G /|p0g| = Tj/,Ij^ u1 g/\Pug\i and take the ratio between the <j> and u> channels in Eq. p2|l . we then have 

R _ zi = r J/4>^4>ti-«PKK = |p i| [xi + yi} 2 _ 3 Q , 13 . 

which is contradicted by the experimental result R exp = 0.15 from BES 0,0] and R exp = 0.75 from PDG [2lj | . 

A solution is to keep the flavour-blind assumption, but allow r ^ 0, i.e. to include contributions from Fig. ^c). 
With the experimental data for the /o(1710) production, we can determine the coupling strength r: 

(1 - Xo)x\ + yi - V2xqZi , (14) 



XoVi - V2zi 



where xo = (2i? e a;p|Pwi|/|P0i|) 1 ^ 2 - Note that R exp is determined by the experimental data; r is not a free parameter 
in this calculation. We find that for R exp = 0.15 ~ 0.75, r has a range of 0.5 ~ 2.5, which suggests a rather strong 



contribution from the doubly OZI disconnected diagram. It shows that though there are significant descrepancies 
between the PDG estimates and BES measurements, they both favor large values for r and breaking of the OZI- 
rule. In particular, this range is within the expectation of non-perturbative glueball-QQ couplings as found in the 
configuration mixing scheme for the fa states, and as suggested by the parameters r-i and in Section II. In this 
sense, the value r ~ 1 turns out to be a self-consistent result of this model. 

To test the self-consistency of the calculation and sensitivity of other branching ratios to r, we examine the following 
four cases: i) with the parameters of Fit-I, and R exp — 0.75 from the PDG [2]]], we derive r = 0.5 and the predictions 
are labelled as Fit-I(a); ii) with the parameters of Fit-I, and R exp = 0.15 from the BES [l3L IT^ |. we derive r = 2.5 
and the predictions are labelled as Fit-I(b); hi) with the parameters of Fit-II, and R exp = 0.75, we derive r = 0.6 
and the predictions are labelled as Fit-II(a); iv) with the parameters of Fit-II, and R eX p = 0.15, we derive r = 2.2 
and the predictions are labelled as Fit-II (b). Recalling that Fit-II has accommodated the data from BES, we note in 
advance that Fit-II(b) is the only set where the new experimental data from BES have been used throughout. We 
hence would expect that predictions from Fit-II(b) are the self-consistent ones. 

With a fixed R exp = T j/^^f^^KK A" J/iji^uf^uKK > we are a ^ e to normalize the glueball production widths via 
Eqs. 1(12(1 and 1)14(1 . In particular, the spin-averaged invariant amplitudes squared are 

\ M <j>G\ - dA^G| - — r- — j— - 1 ; ; —7= - , (15) 

where the ratio br^i^ K ^ = r(/o(1710) — > K K)/T( fo(1710) — > all) will be estimated based on the present experi- 
mental information. So far, the experimental evaluations of br f i^ K g- still have large uncertainties. Earlier analyses 
from Ref. [2J| gave brfi^ K j^ = 0.38, which however may not be reliable. Taking into account that KK is the 
dominant decay mode for the fn(1710), and the branching ratio fractions of other pseudoscalar pair (irir, 7777, 7777') 
decays over KK are small 0,03' a more reliable estimate could be brf\^ K ^ ~ r(/ (1710) — > K K)/T(f (1710) — > 

(KK+Tnr+rjrj+rjr)')) ~ 0.60 with the WA102 data. We have also assumed that the 4-7T decay channel is negligible |2^ . 

A large value of br t\^ K j^ — 0.60 will predict different glueball production rates. However, we would like to point 
out in advance that the magnitudes of the predicted glueball production rates and the predicted branching ratios for 
the /o(1500) and /o(1370) will not change drastically with the change of brfi^ K ^ from 0.38 to 0.60, or even up to 
0.8. This is natural since the underlying physics in the exclusive process of J/ip — > V/q should not be sensitive to 
the experimental uncertainties arising from the total widths of the /q. We will come back to this point in the next 
Subsection. To proceed, we will then adopt brji_ >K g : = 0.60 in the calculations. 



C. Predictions for the glueball, /o(1500) and /o(1370) 



With brti^KK = 0.60, we have access to the normalized invariant amplitudes squared for the glueball production 
via Eq. 1(15(1 . The branching ratios for J/ip — > VG can then be derived, and the results are presented in Table ITVI for 
the four circumstances defined above. 

Although there exist uncertainties due to dependence on the experimental data for br^i^ K ^, we find that the 

magnitudes are well constrained at order 10~ 4 . Interestingly, the partial width of J/tp — > ujG is predicted to be 
larger than J/ip — > <\>G by about a factor of two. A well-determined value for R exp will improve the accuracy of 
this prediction. Meanwhile, it is interesting to compare this approach with the one proposed in Ref. ^^|. In the 
latter, the pQCD assumption provided a way to connect the electromagnetic decay processes with the strong decays 
via gluon-counting rules to the next-to-leading order, through which the glueball production rate can be estimated. 
However, the estimate strongly depended on the validity of treating all the exchanged gluons as perturbative. In 
contrast, the validity of pQCD treatment in this work can be examined by the factorization scheme and controlled 
by the experimental data. It is worth noting, since the doubly OZI disconnected processes are intimately related to 
the glucball-QQ mixings, it is not surprising that non-perturbative contributions are still significant in J/ip — + V/q. 

As follows, we will examine the predictions from those four circumstances for the production of the /o(1500) and 
/o(1370), and compare them with the available experimental data from BES. 

To calculate the partial decay widths for J/tp — > V/q — ■> VPP, we still need the branching ratios for /g — > PP. 
For /q — » KK and irn, the experimental data from PDG are listed in Table Ivl Some of these values still have large 
uncertainties. However, we note that the branching ratio fractions, br(J/ip — * cfyj'^ — + (j)PP)/br(J/if) — > w/q — * loPP), 
are independent of the branching ratios for /q — * PP. Therefore, the theoretical predictions for the J/ip decay 
branching ratio fractions can still be compared with the experimental data. Again, we note that the predicted physics 
for the exclusive J/ip — > V/q should not be dramatically sensitive to the uncertainties arising from the branching 
ratios for /q — > PP. 



As summarized earlier, the BES results [l3| show that the branching ratios for J/tp — > 4>fo ~ > 4> n7r are sizeable 
for /o(1370) and /o(1500), while the /o(1710) is negligible. In J/tp — > (pffc — > <pKK, the situation changes, and the 
/o(1710) has the clearest signals. On the other hand, signals for the /o(1370) and /o(1500) in witty and ujKK are 
found to be small and difficult to isolate, while they are relatively clearer in (pinT and <pKK . 

In Table IVll the branching ratios for the J/ip decays into <pKK, ujKK, (firm, and ujtttt via the /q are presented 
for Fit-I(a), Fit-I(b), Fit-II(a), and Fit-II(b). It shows that Fit-I(a) and Fit-II(a) with R exp = 0.75 lead to sizeable 
branching ratios for the /b(1500) in ujtttt. Therefore, they are not favoured by the BES observations. This should not 
be surprising since BES is in favor of R exp = 0.15. 

In contrast, Fit-I(b) and Fit-II(b) with R exp = 0.15 produce consistent patterns for /o(1500) and /o(1370). The 
/o(1710) branching ratios in both (uj)tttt and (4>)tttt are found significantly smaller than in (uj)KK and (<j))KK, 
respectively. With a large branching ratio of (13.2 ± 2.6) x 10 -4 for J/ip — > o;/o(1710) — > uoKK as input, we obtain 
the corresponding ratio 1.45 x 10~ 4 for J/ip — > w/o(1710) — * ujtttt. This quantity seems quite big, and implies that 
it could be relatively easier to separate the /o(1710) signal in the ujtttt channel than both /o(1500) and /o(1370). 
However, taking into account that the tttt/KK ratio still has large uncertainties, e.g., BES gives an upper limit of 
0.11 with 90% CL, the branching ratio of J/ip — > w/o(1710) — > ujtttt can be further suppressed with smaller tttt/KK 
ratios. The Jo (1500) has sizeable branching ratios in ((P)tttt, while all the other channels are small, especially the ones 
in association with uj production. It is worth noting that an upper limit on J ftp — > w fo (1500) is given by BES at 
0.9 x 10~ 4 [24]]. Adding the branching ratios of (uj)KK and (uj)tttt for the /o(1500), the ratios are well below the 
limit, which is in agreement with the data. Comparing Fit-I(b) with Fit-II(b) for the /o(1370), we find that the large 
uncertainties arising from the KK to tttt ratios do not change the major charactor of the branching ratio pattern. Its 
branching ratios in ((P)tttt are significantly bigger than in (uj)tttt, which was previously regarded as a "puzzle" for this 
state. 

As manifested by Eq. (|14f> . a small value for R exp such as 0.15, corresponds to r ~ 2.2, which implies strong non- 
perturbative contributions from the doubly disconnected processes. We stress again that r is not a free parameter in 
this model. It is determined by the experimental measurement of R exp for the /o(1710), while all the other variables 
in Eq. (|14|1 are from either kinematics or the independently-determined mixing matrix elements. Therefore, the ratio 
R exp also provides a good test of r, which is in a reasonable range between 0.5 to 2.2 in line with strong QCD 
( "non-perturbative QCD" ) . 

Interestingly, we also find that for much smaller R exp , hence larger r, the production branching ratios for /o(1500) 
and /o(1370) via ujtttt will increase and worsen the predictions. In this sense, we identify that the strong contributions 
from the doubly OZI disconnected processes are not an artificial effect. The preferred coupling strength is consistent 
with the expectation of a nonperturbative glueball-QQ interaction appearing in the configuration mixings. It also 
suggests that an accurate measurement of R exp for the /o(1710) [l3L 111] will be crucial for the study of the reaction 
mechanisms. 



IV. SUMMARY AND PERSPECTIVE 

In this work, we proposed a factorization scheme for the study of the scalar mesons /q production in the J/ip 
hadronic decays into the isoscalar vector mesons V (— <p an d uS) and pseudoscalar meson pairs PP (= tttt, KK, rjr), 
and r/rj'), i.e. J/ip — > Vffa — > VPP. With an improved treatment for the glueball-QQ mixing in the configuration of 
the /o, we found that the reaction mechanisms were very sensitive to the structure of the /q. In particular, due to 
the glueball-QQ mixing, the predicted branching ratios for J/ip — ► VJq — ► VPP exhibit unusual patterns, which are 
in agreement with the experimental data from BES. 

The importance of glueball-QQ mixing is also highlighted by the indispensible contributions from the doubly 
disconnected processes, which turn out to be nonperturbative and violate the OZI rule. Since the coupling gg — ► QQ 
in the doubly disconnected processes is essentially the same as the glueball-QQ mixing, the nonperturbative feature 
of the doubly disconnected processes is self-consistent with the proposed configuration mixing scheme for these three 
fo states. In this sense, our results not only provide an understanding of the recent "puzzling" experimental data 
from BES [lii. IT4I llfij. but also highlight the strong possibility of the existence of glueball contents in the /o(15QO), 
and its sizeable interferences in /o(1710). Furthermore, due to the configuration mixing, the \nn) dominant fn( 1370) 
tends to have a lower mass lower than 1370 MeV, which also agrees with a recent more refined analysis MM 

With Eqs. © and (|10|> . and applying the method of Ref. @, the the relative decay widths (excluding phase space) 
for /£->■ 77 are found to be / (1370) : / (1500) : /q(1710) - 12 : 2 : 1. The results are consistent with those of 
Ref. [fj, which should not be surprizing since the mixing matrices are similar to each other. 

In this factorization scheme, a quantitative normalization of the scalar glueball production rate in the J/ip — > VG 
is also accessible. With the pure glueball mass in a range of 1.46 - 1.52 GeV, we obtain the branching ratios 
brj/^p^^Q ~ \br j i^^^g ~ (1 ~ 2) x 10~ 4 . Although a direct measurement of the glueball production seems 



impossible, the success of this approach provides a key for understanding the recently claimed signals for another 
scalar /o(1790). This state is found distinguished from the /o(1710) based on two observations [THl 12^ : i) the 
/o(1710) has clear signals in (uj)KK and (<fi)KK, while the /o(1790) is not visible; ii) there is a definite peak in (jmir 
for the /o(1790), while the signals for /o(1710) is negligible. As a result, the branching ratios between KK and mr 
are larger than a factor of 20. These distinguished features suggest that the /o(1790) may be a radial excited state 
dominated by \nn) which a ppe ars not to be strongly mixed with the ground state. 

As pointed out in Ref. |23, the decays J/ip — ► 7/0 — * 77U, where V = cj, <f>, and p°, can provide further 
independent information about the f states. Using the discussions of Ref. @, we find the widths for J/tp — > 7/3 
satisfy /o(1710) > /o(1500) > /o(1370), which is also consistent with the BES data. Further improved statistics at 
BES for these fo states in both hadronic and radiative decays should be able to quantitize the model predictions, and 
provide us with more information about their structures. 
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7 2 (/o -» rm) [a 2 z l + V2f3 2 y, + r 2 x, + {s/2a - P) 2 r 3 x t ] 2 

7 2 (/o W) 2[a/3(* - v^Wi) + (V2a - /3)(>/2/9 + a)r 3 a:<] 2 

7 2 (/o^^) 3[ Zi + r2Xi] 2 

7 2 (/^ -> ifK) 4[§(* + v/2yi) + 

TABLE I: The theoretical reduced partial widths. Parameters a and /? are defined as a = (cos0 — \/2sin0)/\/3 and /3 = 
(sin <j> + \pl cos (f>)/y/3, where <f> is the flavour octet-singlet mixing angle between r/ and r/' . Factor j 2 denotes the invariant decay 
couplings for /o — * PP. We adopt the same form factors as Ref. 0,0 in the calculations. 



B.R. fractions 


Data from WA102 


Fit-I 


Data of WA102 plus BES 


Fit-II 


D 


2.17±0.9 


2.34 






1/D 






0.08 ± 0.08 


0.07 


/o(1370)^>7>7 
/o(1370)-»KK 


0.35 ± 0.21 


0.62 


0.35 ± 0.21 


0.17 


/ (1500)^7r7r 
/ (1500)^^r, 


o.o ± U.84 


4.09 


0.0 ± U.84 


o no 


f (1600)->KR 
/ (1500)-n™ 


0.32 ±0.07 


0.41 


0.32 ±0.07 
0.78 ±0.66 


0.42 


/ (1500)-» W ' 
/o(1500)^T7J7 


0.52 ±0.16 


0.68 


0.52 ±0.16 


0.73 


/q (1710) -mitt 
/ (1710)-» K 


0.20 ± 0.03 


0.19 


< 0.11 


0.13 


/ (1710)-n,^ 
/ (inO)-fJC/f 


0.48 ±0.14 


0.22 


0.48 ±0.14 


0.20 


/ (1710)-»»,»j' 
/ (1710)— T7J7 


< 0.05 (90%cZ) 


0.03 


< 0.05 (90%cZ) 


0.03 


x 2 




4.6 




3.6 



TABLE II: Fitting results in comparison with the data from WA102 and BES. Quantity D denotes f °^°)^KK ■ The differences 
between these two experimental sets applied for the fittings are highlighted by the bold font. 




Parameters 


Fit-I 


Fit-II 


4> (degree) 


-18.5 ±3.1 


-22.7 ± 3.5 


M G (MeV) 


1464 ± 47 


1519 ±41 


M S s (MeV) 


1674 ± 10 


1682 ± 4 


M nfi (MeV) 


1357 ±23 


1304 ± 7 


/ (MeV) 


84 ± 14 


71 ±5 




1.10 ±0.02 


0.88 ±0.12 




1.12 ±0.27 


1.08 ±0.20 


M 3 (MeV) 


1275 ± 34 


1260 ± 13 



TABLE III: Parameters for the two sets of fittings. 



Br. for G (xlO -4 ) 


Fit-I(a) 


Fit-I(b) 


Fit-II(a) 


Fit-II(b) 




3.97 


1.36 


3.96 


1.63 




9.02 


3.09 


9.11 


3.75 



TABLE IV: Normalized branching ratios for the glueball production in J/tp — > <f>G and J/tp — » u>G. The four sets of calculations 
correspond to Fit-I(a): with the parameters fitted by Fit-I, and R exp = 0.75 from the PDG 21]; Fit-I(b): with the parameters 
fitted by Fit-I, and R exp = 0.15 from the BES 013; Fit-II(a): with the parameters fitted by Fit-II, and R exp = 0.75; and 
Fit-II(b): with the parameters fitted by Fit-II, and R exp ~ 0.15. We note that Fit-II(b) is the one treating the experimental 
data self-consistently. 





KK 


7T7T 


/o(1710) 


0.60 


0.60 x 11% 


/o(1500) 


0.086 


0.349 


/o(1370) 


0.20 x 10% 


0.20 



TABLE V: Experimental values for the branching ratios of /o — > PP. For 6rj ( 1710 )^ 7r ^ and br^^^^K g , we apply the 
results from the BES Collaboration [ Hll4| . 



Fit-I(a) 


KK (<f>) 


Data 


KK (lo) 


Data 


7T7T ((f)) 


Data 


7T7T (lo) 


/o(mo) (xio - 




3.60 


[3.6± 0.6] 


4.80 


[4.8± 1.1] 


0.40 


- 


0.53 - 


/o(1500) (xlO" 




0.19 


(0.8± 0.5) 


1.38 


- 


0.77 


(1.7± 0.8) 


5.60 - 


/o(1370) (xlO" 


- 4 ) 


0.00 


(0.3± 0.3) 


0.08 


- 


0.01 


(4.3± 1.1) 


0.83 - 


Fit-I(b) 


KK (cj>) 


Data 


KK (lo) 


Data 


7T7T ((f)) 


Data 


7T7T (lo) 


/o(1710) (xlO~ 




2.00 


(2.0± 0.7) 


13.20 


(13.2± 2.6) 


0.22 


- 


1.45 - 


/o(1500) (xlO~ 




0.66 


(0.8± 0.5) 


0.09 


- 


2.67 


(1.7± 0.8) 


0.35 - 


/o(1370) (xlO~ 




0.19 


(0.3± 0.3) 


0.05 


- 


1.93 


(4.3± 1.1) 


0.54 - 


Fit-II(a) 


KK (<j>) 


Data 


KK (lo) 


Data 


7T7T ((f)) 


Data 


7T7T (lo) 


/o(1710) (xlO~ 




3.60 


[3.6± 0.6] 


4.80 


[4.8± 1.1] 


0.40 


- 


0.53 - 


/o(1500) (xlO~ 




0.23 


(0.8± 0.5) 


1.20 


- 


0.92 


(1.7± 0.8) 


4.85 - 


/o(1370) (XlO" 


- 4 ) 


0.01 


(0.3± 0.3) 


0.18 


- 


0.14 


(4.3± 1.1) 


1.81 - 


Fit-II(b) 


KK (<p) 


Data 


KK (lo) 


Data 


7T7T ((f)) 


Data 


7T7T (lo) 


/o(1710) (xlO~ 




2.00 


(2.0± 0.7) 


13.20 


(13.2± 2.6) 


0.22 


- 


1.45 - 


/o(1500) (xlO~ 




0.46 


(0.8± 0.5) 


0.13 


- 


1.89 


(1.7± 0.8) 


0.52 - 


/o(1370) (xlO" 




0.26 


(0.3± 0.3) 


0.09 


- 


2.63 


(4.3± 1.1) 


0.94 - 



TABLE VI: Predictions for br 7 ^_, v p—,v pp in comparison with the available experimental data. The BES results [H^.ll4| are 
quoted in the round brackets, while the PDG results |21| are quoted in the square brackets. The symbol '-' means that signals 
of the corresponding states have not been observed in BES experiment. See the caption of Table llVl for definition of these four 
sets. 
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FIG. 2: Process involves strong mixing of nn-ss in the scalar production. This may occur at similar strength to the singly 
disconnected Fig. Q^b) . 



